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SUMMARY  ABSTRACT 


LIDAR  (light  detection  and  ranging.  Laser  Radar)  is  a  transceiver 
system  which  measures  the  amount  of  light  returned  to  a  receiver  due  to 
backscattering  of  the  transmitted  signal  from  aerosols  blocking  the 
path  of  propagation.  The  term  aerosol  refers  to  any  physical  matter 
suspended  in  the  atmosphere,  other  than  standard  atmospheric  gases  (0^, 
^  CC^,  etc.).  All  such  matter,  whether  solid  or  liquid,  possesses 
the  ability  to  reflect  and/or  absorb  light  energy  which  is  incident 
upon  it. 

The  extreme  importance  of  LIDAR,  as  a  tool  for  meteorological  pre¬ 
diction  and  monitoring,  is  that  light  energy,  returned  (backscattered) 
to  the  LIDAR  receiver,  provides  an  observer  with  a  multitude  of  dis¬ 
tinct  information  about  the  volume  of  atmosphere  which  caused  the  light 
scattering,  LIDAR  also  provides  a  remote  ranging  capability. 

This  paper  will  be  devoted  to  an  explanation  of  how  this  returned 
signal  can  be  filtered  for  extremely  detailed  information  concerning 
the  portion  of  atmosphere  with  which  it  interacted.  The  study  was  con¬ 
ducted  in  three  phases.  The  first  phase  is  devoted  to  the  LIDAR 
instrument  itself,  to  include  the  physical  parameters  and  governing 
equations,  as  well  as  the  LIDAR  system’s  advantages  and  shortcomings, 
in  comparison  to  conventional  radar  systems.  The  second  phase  is 
devoted  solely  to  a  typical  vtilume  of  the  atmosphere,  its  physical 
properties,  and  what  must  be  known  about  the  physical  properties.  The 


third  phase  will  be  devoted  to  linking  the  principles  of  the  LIDAR  to 
the  properties  of  an  irradiated  volume  of  the  atmosphere  in  an  attempt 


to  obtain  a  remote  picture  of  the  atmosphere  through  the  characteris¬ 
tics  of  the  returned  signal.  In  particular,  the  measurement  of  number 
density  of  particles  within  a  given  volume,  or  concentration,  is 
desired. 
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In  the  field  of  meteorology,  the  formalism  for  large-scale  atmos¬ 
pheric  phenomena  is  well  understood  and  documented.  The  worldwide  grid 
of  weather  observation  stations,  as  well  as  satellites,  provides  direct, 
measurement  or  observation  of  the  state  of  the  atmosphere  to  a  high 
degree  of  accuracy  and  reproducibility,  provided  the  sc^ie  is  large 
enough.  It  is  the  elusive  micro  and  mesoscale  phenomena  which  yet 
escape  the  interested  observer.  The  minute  scale  most  assuredly  masks 
the  phenomena  within  the  larger,  more  gross  patterns. 

Thus,  an  indirect  observation  approach  must  be  taken  as  is  done  in 
atomic  and  nuclear  physics.  The  structure  of  the  nucleus  of  an  atom 
cannot  be  observed,  or  its  properties  measured,  by  the  use  of  a  stand¬ 
ard  laboratory  microscope.  The  scale  size  is  too  small.  Therefore,  a 
more  indirect  approach  is  used.  That  is,  a  substance  in  question  is 
bombarded  with  well-known  elements,  such  as  protons  or  electrons.  The 
effect  of  the  substance  in  question,  on  the  known  particles,  is  then 
observed  in  an  attempt  to  make  deductions  concerning  the  properties  or 
structure  of  the  unknown.  The  existence  of  mu,  pi,  and  k"  mesons  is, 
for  example,  a  result  of  this  indirect  type  of  observation.  Similarly, 
the  electromagnetic  energy  transmitted  by  a  LIDAR  system  can  be  used  to 
deduce  indirectly  the  properties  of  a  medium  (the  atmosphere)  through 
which  it  propagates. 


PHASE  1.  THE  LIDAR  INSTRUMENT,  ITS  PROPERTIES  AND  GOVERNING  EQUATIONS 

The  UDAR  principle  is  identical  to  that  of  a  conventional  radar 
system  in  that  a  pulse  of  energy  is  transmitted.  For  LIDAR,  the  trans¬ 
mitted  signal  has  a  wavelength  in  the  visible  or  near  visible  region  of 
the  electromagnetic  spectrum  (0.2  to  1.0  microns,  normally).  Typically, 
an  energy  per  pulse  is  on  the  order  of  a  few  joules.  The  duration  of 
the  pulse  is  usually  between  one  and  ten  microseconds,  although  pulses 
shorter  than  one  microsecond  are  desirable  for  ranging  resolution.  The 
signal,  as  it  propagates  through  the  atmosphere,  will  diverge  somewhat, 
depending  on  the  optical  properties  within  the  laser.  Typically,  LIDAR 
beamwidths  vary  from  1  to  20  milliradians.  The  signal,  as  it  propa¬ 
gates  away  from  the  transmitter,  is  scattered  by  the  medium.  A  small 
fraction  of  the  signal  will  be  scattered  in  the  backward  direction.  It 
is  this  backscattered  energy  (0  =  180°,  Fig.  1)  which  is  detected  and 
measured  as  a  function  of  range,  or  distance,  at  which  the  scattering 
occurred.  Two  quantities  are  now  available  to  the  observer: 

a.  The  amount  of  energy  returned  due  to  scattering  from  a  partic¬ 
ular  volume  of  the  atmosphere. 

b.  The  range  at  which  this  scattering  occurred. 

The  manner  in  which  the  energy  is  detected  and  ranged  is 
straightforward.  A  narrow  column  of  the  atmosphere  is  illuminated,  and 
the  light  is  scattered  in  all  directions,  as  well  as  absorbed  by  the 
medium.  Whatever  fraction  of  the  original  signal  that  is  returned  to 
the  detector  is  passed  through  a  narrow  band  filter  to  reject  unwanted 
stray  light  of  different  wavelengths,  such  as  sunlight.  The  signal  is 
then  focused  internally  on  a  series  of  photomultipliers.  The 
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photomultiplier  tubes  convert  the  focused  light  into  an  electrical 
signal  which  is  then  amplified  to  a  usable  level,  calibrated,  and  por¬ 
trayed  in  some  form  of  scope  display  or  digital  readout.  Frequently, 
the  returned  signal  is  so  weak  that  individual  photons  are  counted. 

The  problem  of  signal  amplitude  measuring  due  to  scattering  from  many 
different  ranges  is  one  of  counting  the  number  of  pulses  from  the 
photomultipliers  in  a  series  of  gated  intervals. 

The  amount  of  power  returned  -  i.e.,  the  fraction  of  the  incident 
power  scattered  in  the  backward  direction  -  is  dependent  upon: 

a.  Hu.  intensity  of  the  incident  signal. 

b.  The  wavelength  of  the  incident  signal. 

c.  The  physical  and  optical  properties  of  the  scattering  media. 

LIDAR,  rather  than  a  conventional  radar,  is  essential  for 

observing  small  particulate  matter  and  gases  as  indicated  by  the  equa¬ 
tion  which  yields  the  efficiency  for  scattering  a  signal.  The  scat¬ 
tering  cross  section  (a)  is  defined  as: 


a 


64nSa6 


o  1 2 

m*~  -  1 
m2  +  2 


CD 


for  scattering  particles  with  radii  sizes  less  than  the  wavelength  of 
the  incident  signal, 
where  a  =  radius  of  scatterer 
A  =  wavelength  of  signal 
m  =  ir  .ex  of  refraction  of  scatterer. 

If  \  =  10  cm  (conventional  radar) 


and  a  =  10  mm  (raindrop) 

6 


then 


lr.-4  2 

10  cm  . 


This  ratio  is  easily  detectable  by  conventional  radar. 

If  X  =  10  cm  (conventional  radar) 

-4 

and  a  =  10  mm  (dust  particle) 

a6  in-34  2 

then  —t  =  10  cm  . 

i4 

This  ratio  would  not  be  detected  by  conventional  radar  since  the  dust 
particle  is  so  small. 


If 

and 

then 


>. 

a 

6 


=  10  cm  (LIDAR) 

-4 

=  10  mm  (dust  particle) 


10‘18cm2. 


The  dust  particle  is  now  detectable  due  to  the  much  shorter  wavelength 
of  the  LIDAR. 

Thus,  a  LIDAR  possesses  the  capability  of  detecting  extremely 
minute  atmospheric  particles,  as  well  as  gases. 


PHASE  2.  PROPERTIES  OF  ATMOSPHERIC  SCATTERERS 


It  was  shown  in  PHASE  1  that,  for  a  given  signal  intensity  (IQ)» 
the  amount  returned  (ID)  due  to  scattering  at  a  range  (R)  would  be  con- 
sidsrably  less  than  Iq.  The  actual  amount  returned  is  dependent  on  the 
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scattering  angle  (0),  the  wavelength  (X),  and  the  properties  of  the 
scatterers.  Therefore,  if  the  transmitted  intensity  (Iq)  is  known,  and 
the  returned  intensity  (IR)  is  known,  some  deduction  concerning  the 
properties  of  the  scatterers  can  be  made.  However,  the  computation  is 
complex. 

The  atmosphere  is  composed  of  many  different  particles  and  gases 
of  many  different  sizes.  The  scattering  cross  section  (efficiency  of  a 
particle  to  scatter  light)  was  theorized  by  Mie  (1908).  His  hypothesis 
was: 

0(0)  =  (ij(0)  +  i2(0)}  (2) 

where  o(0)  is  the  scattering  cross  section  for  a  direction  0;  ij(0)  and 

2n 

i2(0)  are  the  amplitude  scattering  functions  (Appendix  I);  and  k  =  -p 
The  Mie  solution  describes  the  ability  of  a  particle  to  scatter 
incident  light  in  any  direction  as  a  function  of  the  radius,  index  of 
refraction  of  the  scattering  particle,  and  the  wavelength  of  the  inci¬ 
dent  light.  An  indication  of  the  completeness  of  the  solution  is  the 
fact  that  it  contains  the  two  classical  cases  of  scattering: 

a.  When  a  particle's  diameter  is  much  less  than  the  wavelength, 
the  Mie  solution  reduces  to  the  familiar  case  (as  in  conventional 
radars)  of  Rayleigh  scattering,  where  a  particle's  ability  to 
scatter  light  is  inversely  proportional  to  the  fourth  power  of  the 
laser  transmitting  wavelength. 


b.  When  a  particle's  diameter  is  orders  of  magnitude  greater  than 
the  wavelength,  the  Mie  solution  approaches  geometric  scattering, 
where  its  ability  to  scatter  light  is  nearly  independent  of  wave¬ 
length  and  is  proportional  to  the  particle  size  or  geometric  area 
of  the  particle. 

2 

a  =  c  a 
o 

In  these  two  regions,  the  variation  of  the  cross  section  (a)  with 
a  and  A  is  quite  smooth.  However,  if  a  is  nearly  equal  to  A,  j  is  a 
highly  fluctuating  quantity  and  the  generalized  version  of  Mie's  equa¬ 
tion  is  required.  From  electromagnetic  theory,  this  equation  can  be 
derived  rigorously  for  spherical  scatterers  of  arbitrary  size.  A  some¬ 
what  detailed  derivation  is  provided  in  Appendix  I. 

An  important  assumption  was  made  within  Equation  2;  that  is,  that 
the  scattering  particles  are  spheres.  The  Mie  theory,  when  applied  to 
scattering  from  nonspherical  particles,  indicates  that  the  shape  of  the 
particle  is  not  nearly  as  important  as  its  size  and  optical  properties. 

The  most  important  quantity  in  determining  the  cross  section, 
other  than  radius  and  wavelength,  is  the  complex  index  of  refraction. 
This  is  a  quantity,  associated  with  the  particle,  which  defines  its 
ability  to  scatter  light;  i.e.,  it  indicates  how  well  a  given  particle 
will  reflect  and/or  absorb  incident  light. 


Now,  given  an  arbitrary  particle  radius  (a),  its  index  of  refrac¬ 


tion  (ra),  an  incident  signal  intensity  (IQ)>  and  wavelength  (X),  the 
amount  of  power  scattered  in  a  direction  (0)  can  be  computed. 

However,  the  typical  atmosphere  consists  of  scattering  particles 
of  various  sizes  and  indices  of  refraction.  Therefore,  the  total 
scattering  cross  section  for  an  illuminated  volume  of  the  atmosphere 
will  be  the  su,.»  of  the  scattering  cross  sections  for  all  the  different 
types  of  particles  (different  indices  of  refraction)  and  the  sum,  over 
the  size  distribution  (the  number  of  particles  for  each  type).  The 
cross  section  for  all  particles  whose  index  of  refraction  is  m  is: 


a 


m 


n(a) am(a)da 


(4) 


CT  = 


(5) 


PHASE  3.  LINKAGE  OF  LIDAR  PRINCIPLES  TO  SCATTERING  PRINCIPLES 


The  returned  power  is  a  function  of  the  incident  power  and  the 
backscattering  coefficient.  The  LIDAR  actually  measures  the  volumetric 
backs cattering  coefficient  (oT);  i.e.,  the  signal  can  be  processed  to 
yield  o^.  However,  these  values  are  not  useful. 

The  process  of  diffusion  in  the  atmosphere,  as  a  function  of  time 
and  space,  is  represented  as  concentration  M. 


M  =  |n  /  a3pn(a)da  (6) 

al 

The  ranging  and  real  time  capability  of  LIDAR  will  automatically 
provide  a  number  which  is  a  function  of  time  and  space,  but  this  ques¬ 
tionable  number  must  be  edited  in  such  a  way  that  it  is  an  accurate 
measure  of  the  concentration  of  particles  (M)  at  any  particular  time 
and  space  coordinate. 

Let  a  tracer  be  introduced  into  the  atmosphere.  This  tracer  shall 
fluoresce  at  a  designated  wavelength,  one  which  the  LIDAR  receiver  is 
tuned  to  detect.  If  the  properties  of  the  tracer  are  such  that  it 
mixes  with  the  atmosphere  and  remains  in  suspension,  then  the  concen¬ 
tration  of  the  tracer  will  be  the  same  as  that  of  the  included  atmos¬ 
pheric  sample.  Since  the  size  distribution  of  the  tracer  is  known, 
in  Equation  4  is  known.  The  index  of  refraction  of  the  tracer  is 
known;  hence,  Equation  5  is  a  constant. 

For  these  conditions  the  quantity  actually  being  measured  by  the 
LIDAR  is  the  total  number  of  tracer  particles,  n^(a).  Hence,  the  con¬ 
centration  of  the  tracer  cloud  is  known.  If  observation  of  diffusion 
of  the  tracer  were  the  only  requirement,  this  method  would  be 
acceptable. 

However,  disseminating  a  tracer  into  a  cloud  each  time  a  concen¬ 
tration  measurement  is  desired  is  inconvenient.  A  tracer  would  be 
totally  unacceptable  if  random  observation  of  air  pollution  is  desired. 

If  a  tracer  is  not  used,  assumptions  must  be  made  concerning  the 
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size  distribution  of  particles  and  integrated  index  of  refraction.  For 
a  typical  atmospheric  sample,  some  particles  would  be  falling  out; 
hence,  the  size  distribution  is  also  a  function  of  time.  The  effect 
of  this  and  variation  of  index  of  refraction  and  size  distribution  on 
the  concentration  are  currently  being  studied  at  the  National  Center 
for  Atmospheric  Research  (NCAR),  No  conclusive  results  have  yet  been 
determined. 

Two  other  methods  could  possibly  circumvent  this  problem  area. 

One  would  be  the  use  of  Raman  scattering.  This  method  totally  ignores 
size  distribution  and  focuses  on  particle  identity.  The  only  short¬ 
coming  of  this  type  of  system  is  laser  technology  itself.  More  power 
output  than  is  now  available  from  commercial  lasers  would  be  needed  in 
order  to  produce  any  useful  ranging  capability.  The  other  solution 
would  be  to  use  a  standard,  but  wavelength  tuneable,  LIDAR  system. 

This  would  provide  more  information  concerning  the  indices  of  refrac¬ 
tion  of  scattering  particles.  Again,  however,  improved  laser  tech¬ 
nology  (power  output)  will  be  required. 

A  LIDAR  using  a  tracer  seeded  atmosphere  would  be  orders  of 
magnitude  better  than  any  diffusion  monitoring  system  which  now  exists. 
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APPENDIX 

DERIVATION  OF  SCATTERING  EQUATIONS  FROM  MAXWELL'S  EQUATIONS 


There  are  four  fundamental  equations  of  electromagnetism — 
Maxwell's  Equations.  The  equations  are  completely  general  and  apply  to 
all  electromagnetic  phenomena  in  media  which  are  at  rest  with  respect 
to  coordinate  system  used.  They  are  valid  for  nonhomogeneous ,  non¬ 
linear,  and  even  nonisotropic  media. 

These  equations  are: 


VxE  +  ||  =  0 


(7) 


-  5 


(8) 


V-B  =  0 


m-  -*• 

V»D  =  p 


(9) 

(10) 


Some  other  useful  relationships  are: 


D  =  eE 


R  =  UH 


(ID 

(12) 


J 


-* 

oE 


(13) 


0 


(14) 


it  +  7,J 

where  E 

= 

electric  field  intensity 

->• 

D 

= 

electric  displacement 

-*■ 

J 

— 

current  density 

H 

= 

magnetic  field  intensity 

-*■ 

B 

= 

magnetic  induction 

P 

= 

free  charge  density 

e 

= 

dielectric  constant 

V 

= 

magnetic  permeability 

0 

= 

conductivity  of  medium 

t 

3 

time 

From  these  eight  equations  the  motion  of  an  electromagnetic  wave 
through  some  medium  can  be  determined.  The  intensity  of  a  beam  at  a 
particular  set  of  coordinates  (x,y,z)  and  at  a  particular  time  ft)  are 
to  be  determined.  Equations  7  through  8  are  rearranged  and  combined 
along  with  Equations  11  through  14  to  yield  the  wave  equation.  Solu¬ 
tions  to  the  wave  equation  provide  the  magnitude  of  E  and  H  fields  as  a 
function  of  x,  y,  z,  and  t.  However,  solving  Maxwell's  equations  in  an 
x,  y,  z,  and  t  coordinate  system  is  a  very  difficult  and  complex 
procedure.  Therefore,  "Plane  Wave  Solutions"  are  defined  as  those 
solutions  to  Maxwell’s  equations  depending  only  upon  one  space  and 
time  coordinate. 

Let  £  be  a  space  coordinate  such  that 
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E  =  fc(S,t) 


H  »  H(?,t) 


SV  n 


_  / 


Figure  2.  Coordinate  System  Illustration. 


where  i  =  unit  vector  in  x  -  direction 


=  unit  vector  in  y  -  direction 
=  unit  vector  in  z  -  direction 


unit  vector  in  5  -  direction 


The  following  mathematical  operators  reflect  the  change  of  the  £ 
axis  with  respect  to  the  x,  y,  and  z  axes: 


Figure  3.  The  5  Coordinate  System. 


In  Figure  3,  the  point  P  has  two  sets  of  space  coordinates 


n  C  ■*■  8E 
7-E  =  n-^ 


Maxwell's  equations  may  now  be  written  in  a  simpler  and  more 


workable  form. 


•*■  3E  3H 
nx—  +  wtt  =  0 

3 1  at 


-  3H  3E  t 
nxH  '  e3t  “  oE  =  ° 


--  3H  3  ,-►  :►  rt 
n‘3£  =  3^n‘H^  =  0 


Since  p  is  equal  to  0  in  a  conductive  media  -  i.e.,  no  free 


charge  -  then  V-D  is  equal  to  0.  Therefore 


a 

~(n*E)  =  0 

ot,  ■ 


Differentiating  Equation  15  with  respect  to  ? 


2*  2* 
-  34E  a  3  H 

""a?2  VStK 


Crossing  n  into  Equation  19 


lV  -*■  32E 
—<nx  nx — =• 

“I  .  H2 


*|Sxfe 


Differentiating  Equation  16  with  respect  to  t 
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0. 


(21) 
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A  plane  wave  solution  of  Equation  23  is 


* 


ikm£  +  i(ut 
e 


where 


k 


tu  2n 
c  ~  *T 


and 


m 


4ttipg 

U) 


Now  that  a  representation  for  r. hn  ;uopagation  of  an  electromagnetic 
wave  through  a  medium  has  been  obtained,  boundary  conditions  are 
imposed;  that  is,  the  effect  of  a  sharp  boundary  between  homogeneous 
media  is  to  be  considered.  For  the  atmosphere,  these  boundaries  are 
the  gases  and  aerosol  scatterers.  An  expression  for  the  change  in  the 
electromagnetic  wave  (scattering)  caused  by  these  boundaries  can  be 
formulated  by  considering  the  solutions  to  the  wave  equation  with 
specific  boundary  conditions. 

The  "disturbances"  (u,v)  of  a  wave  can  be  expressed  in  terms  of 
the  plane  wave  solutions  of  Equation  23.  This  solution,  as  shown  by 
Mie,  gives  the  components  of  the  E  and  H  fields  at  any  point  inside  and 
outside  a  scattering  particle. 


2n+l 

(n+1) 


Pn(cos©) 


(24a) 
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Figure  4.  Scattering  Diagram. 


Consider  a  simple  example  of  this  scattering  phenomena,  as  illus¬ 
trated  in  Figure  4.  An  isolated  particle  is  irradiated  by  a  plane 
electromagnetic  wave  (as  previously  defined)  of  intensity  Iq.  For 
propagation  of  this  wave  through  a  homogeneous  medium  (e^p^o^),  its 
character  is  preserved;  however,  the  presence  of  a  scattering  particle 
with  properties  differing  from  the  surrounding  medium  (e2,u2,o2)  dis¬ 
torts  the  wave. 

The  disturbance  of  the  incident  wave  is  classified  according  to: 
a.  That  portion  of  the  incident  energy  which  is  scattered  or 
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reflected  by  an  angle  0. 

b.  That  portion  of  the  incident  energy  which  is  absorbed  by  the 
scattering  particle. 

The  scattering  capability  of  an  arbitrary  scattering  particle  is 
defined  by  its  optical  properties;  i.e.,  its  complex  index  of 
refraction, 

y4iriuo 
-  —r- 

where  the  complex  term  defines  the  absorption  properties  of  the 
scatterer. 

From  the  disturbances,  as  defined  by  Maxwell's  equations  (Equa¬ 
tions  24a  and  24b),  a  set  of  usable  amplitude  scattering  functions 
(8^(0),  S2(0),  S3(0),  S4(0))  are  defined. 

00 

SAO)  =  21  -v— ■  iv  {a  it  (cos©)  +  b  r  (cos©)}  (25a) 

1K  J  ^  n(n+l)  nnv  nnv  1  K 


s2(0) 


=  2-  —? - TV  (a  11  (COS0)  +  b  T  (cos©)} 

,  n(n+l)  n  n'-  1  n  n*- 

n=l  v 


(25b) 


For  a  scattering  medium  which  is  spherical  and  homogeneous,  the 
functions  3^(0)  and  S4(0)  are  zero. 


’n(cosS)  =  1IHePntcose) 
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Consider  an  atmospheric  model  consisting  of  one  scattering  sphere 
in  a  medium.  For  incident  radiation  of  intensity  1^,  the  radiation 
scattered  in  a  direction  0  (as  illustrated  in  Figure  4)  is  given  by  Mie 
as: 


1(0) 


ij(0)  +  i2(®) 


.2  2 

k  a 


(26) 
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where 


k  -  — 

k  T 

ij(0)  =  ISj  CO) |2 

i2(0)  =  |S2(0)|2 

Hence,  Equation  26  is  a  general  representation  of  the  intensity  of 
radiation  scattered  into  an  angle  0,  due  to  the  presence  of  a  scat¬ 
tering  medium  as  defined  by  i^(0)  and  i2(0).  The  quantities  ij  (G)  and 
i2(0)  are  called  the  intensity  functions. 

The  efficiency  of  an  object  for  scattering  radiation  back  to  the 
source  is  expressed  as  the  radar  cross  section  (o)  of  the  object. 
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i  ■  Radar  (LIDAR)  is  a  transceiver  system  which  measures  the  amount  of  light 

returned  to  a  receiver  due  to  backscattering  of  the  transmitted  signal  from  the  inter¬ 
vening  media  (gases,  droplets,  and  aerosols).  Besides  providing  a  ranging  capability, 
LIDAR  is  an  important  tool  for  the  monitoring  of  atmospheric  diffusion  processes. 

The  usefulness  of  such  a  system  is  most  evident  in  light  of  the  fact  that  back- 
scattered  energy  which  is  detectable  to  the  receiver  provides  an  observer  with  a  mul¬ 
titude  of  information  about  the  small  volume  of  atmosphere  which  caused  the  scattering. 
Since  several  highly  variable  atmospheric  parameters  define  the  amount  and  manner  of 
light  scattering,  hence  the  amount  of  energy  returned,  numerical  models  of  the  atmos¬ 
phere  are  constructed  which  account  for  the  fluctuations  of  the  returned  signal,  in 
terms  of  the  controlling  atmospheric  parameters.  Thus  the  amplitude  of  the  signal 
returned  is  an  anomalous  measurement  of  the  state  of  the  atmospheric  volume  which 
caused  the  scattering. 

In  order  to  relate  the  atmospheric  variables  to  the  LIDAR  system  variables,  an 
exact  equation  which  describes  the  motion  and  interaction  of  the  signal  beam  with  the 
atmosphere  must  be  developed  and  modeled..  The  development  of  equations  of  motion  which 
describe  the  propagation  of  an  arbitrary  electromagnetic  wave  through  a  media  was  first 
accomplished  by  Mie  in  1908  and  later  by  Van  de  Hulst.  This  derivation,  beginning  with 
Maxwell's  equations,  and  appropriate  solutions  to  the  wave  equations  for  the  purpose  of 
a  radar  system,  is  provided. ^  A  discussion  of  the  implications  and  significance  of  the 
scattering  equations  and  parametric  relationships  is  then-  undertaken,  in  an  attempt  to 
isolate  concentration  of  the  atmospheric  volume  as  a  measurable  quantity  in  terms  of 
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